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Abstract 

Burdzy, Pal, and Swanson [2] considered sohd spheres of small ra- 
dius moving in the unit interval, reflecting instantaneously from each 
other and at x = 0, and killed at rr = 1, with mass being added to the 
system from the left at rate a. By transforming to a system with zero- 
width particles moving as independent Brownian motion, they derived 
a limiting stationary distribution for a particular initial distribution, 
as the width of a particle decreases to zero and the number of particles 
increases to infinity. This space-removing transformation has a direct 
analogy in the isomorphism between a new unbounded-range exclusion 
process and a superimposition of random walks with random bound- 
ary. We derive the hydrodynamic limit for these isomorphic processes, 
demonstrating that this elastic exclusion is an appropriate model for 
the reflecting Brownian spheres in one dimension. 



Preliminaries 

This paper is primarily concerned with the hydrodynamic limit of a exclusion 
process Zt on a bounded one-dimensional lattice of grid size with the 
following dynamics: a particle p moves into an adjacent unoccupied site at 
rate proportional to the size of the block of occupied particles of which p is 
a member. We think of each particle in the block as having internal energy 
transferred to the outermost particle elastically. In addition, the process is 
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boundary driven: at constant rate, the leftmost block of particles (possibly 
empty) is shifted to the right one position, and a new particle is added to 
the vacant first position. Finally, particles are killed when they move to the 
rightmost site. 

We call the model boundary- driven elastic exclusion. We are interested 
in the limiting shape of the empirical distribution for all times as the grid 
size scales to zero and the dynamics scale appropriately, and in theorem 5.1 
we prove that this hydrodynamic limit satisfies the differential equation 



with appropriate boundary conditions. This particle system was chosen to 
approximate the system of one-dimensional crowded Brownian spheres de- 
fined in p], and the connection is of interest because the hydrodynamic limit 
of the exclusion process Zt matches the conjectured hydrodynamic limit of 
the Brownian process in that paper. Because of the connection, and because 
the method of proof in that paper also describes the key isomorphism that 
we use to derive the limiting equation, we briefly describe that process and 
the key transformation. H. Rost [12] also considered reflecting Brownian in- 
tervals and derived the hydrodynamic limit above in the case of the entire 
real line. 

Consider intervals = {B^,B^ + 1/A^), such that when > and 
\B^-Bi I > 1/A^ for all j such that j ^ k, B^ moves as independent Brownian 
motion. Intervals 1^ reflect instantaneously and symmetrically, and are killed 
when B^ + 1/A^ = 1. Finally, for k greater than some fco, Bq = —{k — ko)/N 
and B^ = Bq + at until B^ = 0, so that particles continuously enter the 
interval at rate aN. In order to derive the limiting stationary distribution 
in the case ko = 0, the authors of [2J consider the transformations Tj : 



Tt maps If to a point C^, and simply translates unoccupied space, so the 
move as independent, symmetrically reflecting Brownian motions with drift 
—adt, due to the continually inserted intervals. Furthermore, Tt(l) = St is 
a random boundary that changes proportionally to the number of particles 




(0.1) 



(-00, 1] -> [0, St] with 
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entering or leaving the system. Since the distribution of symmetrically re- 
flecting Brownian motions is identical to that of independent particles, we 
are reduced to the case of independent dA'j: = dW^ — adt, reflecting at and 
killed at St- This leads us to conjecture the following hydrodynamic limit, 
which is a form of the well-known Stefan melting-freezing problem: 

Definition 0.1. Given sq > 0, fo G C^([0, sq]) with sq = 1 — Jq° vo{x)dx, 
and a > 0, a pair {v,s) such that s G C^{[0,T]), s(0) = sq, s > 0, and 
V e C^iDr) n C\D^), where Dt = : < x < s(t),0 < t < T}, 

satisfying 



dtv{x, t) = dxxv{x, t) + adxv{x, t) 


< X < s(t),t > 0, 


(0.2) 


f (x, 0) = fo 


< X < s(0). 


(0.3) 


{dxv{x, t) + av{x, t))\x=o = -a 


t > 0, 


(0.4) 


v{s{t),t) = 


t > 0, 


(0.5) 








s{t) = 1 — f (x, t)dx 
Jo 


t > 0, 


(0.6) 



is called a solution to the Stefan problem \0.Sti0.6[ with initial data {vq, sq). 

Condition (10.61) is more familiar in the differential form s'(t) = —a — 
dxu{x,t)\s(^t)- The Stefan problem has been well studied in many forms, 
though perhaps not this exact form. The book [Tl] by Meirmanov is an ex- 
cellent reference. In particular, with a = 0, this equation is the classical one- 
dimensional, one-phase melting probem, where in the region < x < s{t), v 
represents the temperature of water above freezing, and x > s(t) represents 
a region of ice with temperature 0. Strong existence and uniqueness of this 
case is covered in Cannon's book [4j. There is every reason to believe exis- 
tence holds for the equation above as well, due to the natural physical model 
and intrisic boundedness, but we do not take it up in this paper. We should 
also note that the condition Vq G C^([0, 1]) may not be necessary, but serves 
only to make the definitions simpler. Indeed, the main theorem below holds 
whenever the initial condition is in and the solution satisfies the integral 
form 13.11 The Stefan problem has been studied in a probabilistic context 
as well, as a hydrodynamic limit by Chayes and Swindle [5], Gravner and 
Quastel |S], Landim and Valle [TD], and Bertini et al P]. In [5], [TD], and 
[T], the particle model is simple exclusion, with different particle types rep- 
resenting the liquid and solid regions. Our model is close to that of Gravner 
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and Quastel, who use the Stefan hydro dynamic hmit of a zero-range process 
to prove shape theorems for internal diffusion-hmited aggregation, but our 
proofs are not similar and the application is different. 

We now describe the second discrete process, Y/^ , which is the discrete 
analogue to the distribution of the transformed Brownian motion process. 
For notational simplicity we will omit from the process, but it will always 
be used in the corresponding probability measure P^. Let 



A 



N 



2j + 1 2N -1 



2N 2N 



Our state space for Yt is the subset Qn of N^^ such that for t] E Qn, rj^ 
counts the number of particles at site x for a distribution of particles on 
with the following restriction: there must be M particles, with M < N, and 
the particles may only occupy sites x = {2j + 1)/2N with j < N — M. Let 
Mt be the number of particles at time t, and define a random boundary 5"^ 
with 

At exponential random times with rate A^^ for each direction, particles move 
as independent random walks, reflecting at the leftmost site. If a particle 
hits St at time t, it is killed (removed from the system), and St = St- + 1/N. 
In addition, there is a drift effect, occurring at rate aN, for a > constant, 
where every particle except those at zero shift one site towards the origin, 
an additional particle is added at 1/2N, and St shifts one site left. The only 
state T] E for which this does not happen is ?7i/2Ar = A^ — 1, in which case 
there is no change (think of the generated particle being immediately killed). 

The sum of delta masses of weight Yt/N at each site gives a measure 
fiYt or just /it, depending on context, of mass less than one. The object 
II. is an element of the Skohorod space of right-continuous paths on the 
metric space of positive measures on [0, 1], D{[0, 1], A^). Let be the 
probability measure on right-continuous paths in Q^- that determines the 
process Yt. Let the corresponding probability measure on D([0, 1],A1) be 
Q'^ . The hydrodynamic limit of the is the subject of our first theorem. 
For technical reasons, it is more natural to state the convergence in terms 
of a process Xt, described below in detail, such that p{Xt) = Yt, where 
p{x) = (a; — 1)V0. Q^j- is the state space for Xjy, in one-to-one correspondence 
with fi^r. A precise version of the following theorem is found at 14.11 



4 



Theorem 0.1. Suppose that for each P , converges weakly to a mea- 
sure with fixed density Uq G i^^([0, 1]). Then the empirical measures of X^^ 
converge to the unique solution of a weak version liS. 1\) of the Stefan problem 
l\0.^) - l[07S\) with initial data uq- If a solution for the problem l\0.^) - l[07S\) with 
initial data (fo, sq) = {p{uo),mf{x : uo{x) = 0}) exists, then the empirical 
measures of the process Yt converges to that solution in probability. 

The weak version of the problem is defined in section 3. The meth- 
ods used to derive the hydrodynamic hmit are largely based on those in 
the book [H] by Kipnis and Landim, so we identify our contribution in two 
areas. First, although the Stefan problem has been well studied as a hydro- 
dynamic limit, the exclusion process we describe and the application of the 
free boundary problem to such a process is new. Hydrodynamics of exclusion 
processes is an active field of research, but the most general results are for 
gradient systems with finite-range interactions, as in [Bj. The interactions 
of Xt have unbounded range. Second, the simultaneous drift effect of the 
transformed process is unusual, but required by the isomorphism. The non- 
standard process and the simple setting in the unit interval allows for an 
interesting application of elementary harmonic analysis for the bound 
and the uniqueness proof. 

Our general approach, following [S], is to make the free boundary go 
away by building it into the zero-range dynamics of a process Xf^ as de- 
scribed below. From the other direction, the differential equation transforms 
into a nonlinear integral equation which mirrors the form of the process. 
The proof of Theorem lU.ll is in four steps. In Lemma II. 2[ we prove that the 
Markov process describes a relatively compact sequence of probability mea- 
sures. Lemma [2.11 guarantees that any limit points lie in L^([0,1] x [0,T]) 
almost surely. Lemma 13.11 shows that such limit points must satisfy a weak 
version of equation fl0.2p - fl0.6p . Finally, Lemma [4.11 proves that the solution 
of such an equation is unique. Combining these results, we see that the pro- 
cess converges to a measure which is the delta measure on the solution of 
an integral form of the problem which coincides with the solution to that 
problem when it exists. 

In section 5, we show that the two discrete processes described in this 
section are in fact isomorphic, and use the isomorphism to prove the hydro- 
dynamic limit f lO.ip in Theorem 15.11 
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1 Construction and relative compactness 



We construct a Markov process , henceforth Xt, by defining its infinites- 
imal generator. For > let An = {1/2N,3/2N, ..{2N - l)/2N} and 
Mn = N"^^. Let Ai be the set of finite measures on [0, 1], and we associate 
Tj E M.N with its empirical measure in J^,fir] = S^gAjv Consider the 

following generator on functions / : A^tv — ?■ R: 

Clfir^) =aN{f{a{r^))-f{r^)), 
where p{x) = (x — 1) V 0, 

{rjy — 1 ior y = X and x + i/N G , 
Vy + l foTy = x + t/N, 
rjy otherwise, 

and 

{riy + rjy+i/N for y = 1/2N, 
for y = (2iV - 1)/2N, 

rjy+i/N otherwise. 

The appendix of [9] describes the construction of a Markov process Xt on 
Aijy from such a generator (such that d/dtE [f{Xt) \ X^ = rj] \t=s = Cf^r])), 
the idea being that states are changed at the minimum of exponential ran- 
dom times with rates N'^p{r]x) and aN, to the corresponding state, with the 
minimum itself being an exponential random time, well defined almost surely. 

Let fl'j^ be the set of states r] such that ^^.g^^ rj^ = N, rj^ > for x less 
than some b and rj^ = for x > b. When Xq G fl'^ a.s., the process p{Xt) is 
equal in distribution to the process Y/^ described in the introduction since 
it has the same dynamics. The refiecting random walk effect is due to the 
fact that a pile at site x loses particles at a rate proportional to the height 
p{Xt{x)), as if each particle is moving independently in each direction at rate 
A^^. We now consider the drift and random boundary. Let St = min{x G 
An : Xt{x) = 0}. When a particle moves from 5*^ — to 5*^, which can only 
happen if Xt{St — > 2, it is killed, in the sense that p{Xt) no longer 
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counts it, and the boundary St is incremented by l/N'. as desired. With one 
exception, when the process moves from r] to (j{ri), p{a{ri^)) = except 
at a; = 0, where p{a{ri)o) = pirjijM) + pijlo) + 1, representing the generated 
particle. The only exception is the state Xt{l/2N) = N, St = 3/2N, in 
which case a{-) has no effect, again as desired. Thus p{Xt) with boundary 
St is identical to Yt in its dynamics, and therefore in distribution, given 
corresponding initial distributions. Finally, note that Q'^ is closed under the 
process, and can function as the state space, corresponding to the state space 

of Yt. 

Next, we calculate the generator applied to a linear functional. First, we 
have 



where and are operators with N x N matrices: 



/-I 1 ...\ 
1-2 1 .. 



V- 



1-21 

1 -ly 



and 



D 



N 



N 



/ 

-110 
0-110 



A 



y... -1 ly 

Here A* denotes the transpose of A. Note that and Dn represent the 
second symmetric and first left difference quotients, respectively, for functions 
/ e C^{[0, 1]) with /'(O) = /'(I) = 0, in that for these functions, D^fix) 
converges to f"{x) as N goes to infinity, uniformly in A]\f. 

We briefly digress to discuss the choice of our class of functions. Through- 
out the paper, the functions / e C^{[0,1]) with /'(O) = /'(I) = will be 
used for test functions, as they serve several purposes. First, as noted, it is 
these functions for which the operators above converge uniformly to /" and 
/', respectively. Second, they are as a dense class of functions in C([0, 1]) 
and can be used to define a metric on A4. Third, it is against these test 
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functions that the weak form of the Stefan problem holds. Finally, the sub- 
family {^/2 cos(-Kkx)} is an orthonormal basis for the discrete and continuous 
domains, and is used to prove essential bounds later in the paper. 
Returning to our calculation, for / : [0, 1] — )■ R, let 

and by linearity oi jCn, 

^{AMf,p{v.))N-a{D^f,rj)N. (l-l) 
Next we prove relative compactness oi X in D([0, T], Al). Let 

{f,fJ') = J f{x)dn{x). 
We define a metric on A4, the space of positive measures on [0, 1], letting 

= f MiiizJiii^, (1.2) 

i=o 

where /-,• are in C^([0,1]) with /'(O) = /'(I) = 0, a set which is dense in 
C([0,1]). When fi{dx) = u{x)dx, we will use {f,u) and (/, /i) interchange- 
ably. Since each fj is bounded, A & J\4 is precompact if and only if 
is bounded over A (each {fj,fJ^) is bounded and converges along a subse- 
quence, which implies subsequential convergence in the metric, and A4 is 
complete with respect to rf). Note that the sTipports of all Qn are contained 
in Ml = {fi : = 1}, a compact set. Convergence is weak conver- 

gence (convergence of expectations of continuous functions) in the space of 
probabihty measures on the Skohorod space D{[0,T], Ai) oi right-continuous 
functions on A4. For our purposes, convergence in this space can be conver- 
gence in Ai, uniformly in t, since our limit points are continuous. Thus if 
/(/i) is continuous on JH, then f{ij,t)dt and supQ<j<2. /(/it) are continuous 
on D{[0,T], M.). By {Xt}J^Q, we will mean the Markov process on A^ with 
probability measures P^. By /if we will mean the corresponding coordinate 
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process on Aii, with probability measure , so that, for example, for A 
Borel, 

P''[{f,X,)^eA] = Q''[{f,f,,)eA]. 

Relative compactness in this space follows from the following conditions, 
found in Chapter 2 of [9] . Let 7t be the space of stopping times of the usual 
filtration, bounded by T. 

Lemma 1.1. Let Q'^ be a sequence of probability measures on D{[0,T], Ai) . 
The sequence is relatively compact (in the sense of weak convergence) if: 

1. For every t in [0,T] and every e > 0, there is a compact K{t,e) C Ai 
such that supj^ Q'^lfit 4- — ^■ 

2. 

limlimsup sup P^[p{p^Ti P^{T+e)hT) > e] = 0. 

We prove the following lemma by checking these conditions. 

Lemma 1.2. For an initial distribution such that G Q'j^ a.s., the sequence 
{Q^} is relatively compact in D{[0,T], Ai) . 

Proof. Note that condition (1) is automatically satisfied since, for all N, 
P'^lfJ't ^ A^i] = 0. To check (2), we determine the square variation process 
for the P^-martingale Mt = {f,Xt)N - CN{f, Xs)Nds for / e C\[Q,1]) 
with /'(O) = /'(I) = 0. Exactly as in the proof of Theorem 3.2 of |3], we can 
show that — Bgds is a martingale, where, on {Xt = rj}, 

Bt = lim (l/s)E^[((/,Xi+,)jv - {f,Xt)l I Xt] 

s— >0+ 

+ aN{{f,a{r,))j,-{f,r^)Nf 
= N' J2 p{v.)^,{DNf{x + ir + D^f{xr) + aj^{D^f,r^)%. (1.3) 

Since \Dr,fix)\ < \\f'\U and 1/iV E.^a, p(r/.) < l/N^.^A^V. = 1, 

\Bt\<UfT + ^\\fT- (1-4) 
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Fix T e Tti and by r + ^ we will mean (r + ^) A T, and 

fT+6» 



AT 



^0 



BAs I Jv 



f 

Jo 



Bsds, 



so 



Now 



B.ds 



< 



ce 



E[{Mr+e-Mr) 



T+0 



jCN{f,Xs)ds 



P[\M^+e - >e]< 



E[Ml 



< 



ce 

TV?' 



(1.5) 



and \ j^^^ CN{f,Xs)ds\ < C6, since the generator is the inner product of 
derivatives of / with measures of bounded mass. Thus, 

P''[\{f,X,^0)-{f,X,)\>e]<C,e. 

To bound the metric by a given e, we only need consider finitely many 
and choose appropriately for each of these. The bound is independent of 
T, so 

sup P^[p{iir, fJ'T+e) > e] < C'eT (1-6) 

TeTr,e<7 

and (2) is satisfied. Thus Q'^ is relatively compact and has subsequential 
limits. □ 



2 Limit measures are almost surely 

In this section, we prove that subsequential limits of the measures are 
uniformly bounded in L^, depending on the norm of the limiting initial 
distribution. This allows us to apply the convergence and uniqueness results 
of later sections. 
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Lemma 2.1. If for each N, under P is a random variable with values 
a.s. in n'j^ such that sup^ -E'^[l/A^ ^^.^^^ Xo(a;)^] < oo, a subsequential 
limit of the corresponding on D{[0,T], Ai) has the property that /i 
is absolutely continuous with density u G L^([0,T] x [0, 1]), Q°°-a.s. 

Proof. We prove Lemma 12.11 by looking at the evolution of a variant of the 
norm. Let ipo = 1 and ijjk = \/2cos{kTcx), and let 0o = and (pk = 



V2sm{k7i{x - 1/2N)). Recall that An = {1/2N, (2j + 1)/2A^, ...{2N - 
1)/2A^} and let il)j^,<p^ represent the vectors of the respective functions 
evaluated on A^- Then {ipk}'kLo has the following properties: 

1. AniI^^ = -4N'^sm\nk/2N)il)^. 

2. DniI^^ = -2Nsm{TTk/2N)cf)^. 

3- {'4'k^}k=o orthonormal basis for R"^^. 

4. {r]i,V2)N = Y.o~^('^k,Vi)N{^k,V2)N for r]i,r]2 E N^^. 

5. = Eo''{^k,v)Nil^j^ for r] G N^^. 



The first three can be easily checked by calculations, and the last two are 
consequences of (3). Let A^^iv = 2A^ sin(7rA;/2A^). Next we consider, for 



We restrict to G Q'^ so that {4'o,ri)N- = 1. For 1 < < A^, we have 
Afc.n > C > 0, independent of k and A^, so 



7] G N 





Apply the generator 
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where the last bound follows from the inequalities N(ijjk(x — 1/N) — ijj]^(x)) < 
Xk,N and {l,p{rj.))N < 1- Also, 

^Ni'^k, r])N = -aN({'ijjk, (t{v))n - (V'fc, v)n) 

< 2aAfe,jv('0fe, 'n)N{(t>k, v)n + >^1,n/N: 

and together we get 

N-l 



k=l 

N-l 



Next, let bk — {ipk,v)N, and consider 



^'''^ k=i 

^ Xk,N 
k=l j=0 ' 



Wc claim that 

M-l 



El 71 



1=0 



tt/cA / 7rj(2i + 1 



cos 



N \ 2N 



1 /7r(j + A;)\ / f27rj + 7rk\ f 27: jM + 27rkM - 7rk\\ 

CSC I I I COS I ~ ~ I COS I ~ I I 

A \ 2N J \ \ 2N J V 2Ar // 

1 fTr{j-k)\ f f 2TTjM -2TTkM + 7rk\ f 2Trj - Trk\\ 

+ 4'"' 27V J 1^'°% ; -"^^^ 27V J J ' 

(2.1) 

and we prove by induction in M. Since the constant terms of the right hand 
side are the variable terms evaluated at M = 1, we see that we can check the 
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difference S{M + 1) — S{M) to obtain a telescoping sum on the right hand 
side. In other words, we require 

. fnkM\ f7ij(2M + l)\ 1 fTT(j + k) 

sm cos — ^ = - CSC ^ 

\ N J V 2Ar J A \ 2N 

f f27rjM + 27rkM -7rk\ f 27: j{M + 1) + 27rk{M + 1) - 7rk\\ 

X I COS I I COS I I I 

\ \ 2N J \ 2N J) 



1 f7r{j-k) 

+ - CSC ' 



4 V 2iV 

/27rj(M + 1) - 27rk{M + 1) + tt/cX /27rjM - 2TTkM + tt/c 



X cos 



I ^ 



2N 



Recalhng that cos{A + B) -cos{A) = -2sin(S/2) sm{A + B/2), first with 
A = {27rjM + 2nkM - 7rk)/2N and B = n{j + k)/N, the first term on the 
right hand side becomes 

1 . f 27rjM + 27rkM + 7rj 



- sm 



/ 27rjM + 2TrkM + ttj \ 



and with A = (27: jM - 27rkM + 7rk)/2N and B = 7r{j - k)/N, the second 
term becomes 

1 . f 27ijM -27ikM + 7ij' 



— sm 
2 



/27rjM-27rA;M + 7rA 
V 27V ) 



Finally, the difference formula for sin is sin(/l+i?)— sin(74) = 2 sin(i?/2) cos(74+ 
B/2), and substituting A = {27rjM - 27rkM + 7rj)/2N and B = 27rkM/N, 
we get the desired formula. 

Therefore, for j — k even, we get 

, , , , 1 fTr{j + k)\f f27rj + 7rk\ f -Tik 



AN \ 2N J \ \ 2N J \2N 
1 fTr(j-k)\ f f Tik\ /27rj-7r/e 

+ TTT CSC — COS — — — COS 



AN \ 2N J \ \2N J \ 2N 

1 f^{3'^k)\ ( / Tik 7ij\ ( i\k 

CSC — — COS — — + — — cos ' 



AN \ 2N ) \ \2N N ) \2N , 

1 fT^{j-k)\f f-Tik\ f-irk 
+ 4]V^^H^i^Ji'°Hw N 
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and using the identity cos(A + B) — cos{A) — —2 sm{A + B/2) sin(S/2), we 
get 

{<Pk,ipj)N = - 7^ (sin {^] + sin ( J = 0. 



2N V K'^N J V . . 

For j — k odd, wc get 

1 /7r(i + fc)\ / f27rj+7rk\ f -irk 

(A, ^,)n = ^ CSC {-^^ j (^cos [-^^ j + cos 

1 f T^ij ~ k)\ f f nk\ ( 271 j — nk 

+ TTT CSC — — COS — — — COS ' 



4N \ 2N J \ \2N J V 2iV 

and this time, with cos(>l + B) + cos(A) = 2 cos(A + B/2) cos(S/2), 

(Mn = ^cos (cot (^^^1^) - cot 

Then substitute 

cot(>l + B)- cot{A) = -2csc(A) sin(S/2) cos(S/2) csc(A + B), 
giving 

{(pk,i'j)N _ COs(^)cOs(f|) 



2A^sin im 27V2 sin (%^) sin (^) ' 
We can conclude that 

bkbj{(j)k,'ipj)N _ bjbk{(j)j,ipk)N 
2iVsin(^) " 2iVsin(|^) ' 

and, after cancehng pairs for j > 0, and recaUing that bo — 1, the double 
sum reduces to 

N-l N-1 , , / \ N-1 , , > 

V— V— OfcOj(PA-- Vj/N _ \ - bkjPk- I:\)}n 

r. Afe AT ~f Afe AT 

fe=l j=0 ' k=\ ' 



N-1 



bka{k) cos (^) 



.=1 2Ar2sin2(^) ' 
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where a{k) = k mod 2. We have 2A^2gjn2 ('^^ > CA;^ for 1 < A; < - 1, 



so 



N-l N-l 

EE 

k=l j=0 



bkbj{(l)k,ipj)N 



fe=l 



Finally, noting that {ri,ri)]\f < {pir].) , r}) n + 1, we obtain an estimate for the 
generator applied to hj^, 

The process hN{Xt) — CMhN{Xs)ds is a martingale, so 

E"" [hN{XT) + 2 j\xt,Xt)Ndt^ < E'^ihNiXo)) + CT, 
and since > 0, and HnIXq) < {Xq,Xq)n, 



E 



N 



{Xt,Xt)Ndt 



<£;^ [{Xo,Xo) n]+CT, 



where C does not depend on N. The proposition below finishes the proof of 
the lemma. □ 

Proposition 2.1. Let — )■ Q'^ he a weakly convergent sequence of proba- 
bility measures on -D([0, T], A^) representing the empirical measures of Markov 
processes {Xt,P'^) onN^^. Suppose that 



sup E 

N 



< oo. 



Then ii{dx,t) is absolutely continuous with density u e L^([0,1] x [0,T]) 
Q'^-a.s. 

Proof. We consider the mollification K^p, defined by 



Ken{x, t) 



2e 



lji*{dx,t), 



where is the projection of n onto the torus T created by identifying and 1. 
We will prove that /i* has density, which implies that /i does. Henceforth 
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we use the notation /i for simplicity. For given N, let e — + 'm/N, where 
< eo < l/AT. Then let 



1 

Ye 



X(-m/N+eo,m/N+eo]d^l{y). 



If fj, is the empirical measure for r] e N^, then we can calculate K^n{x)'^dx. 



Kel^{x) 



2eN 



2m 



x—eo,x-eo+ 



—m . 

N 



xeAf 



k=l 



So we calculate 



1 Z^"" 



< 



1 

1 

2 

m 



,jfc=i 



2m- 1 



x&Ai, 



x&Ai^ k=l 
' 2m-l \ 



-E 

N ^ 



Vx I ' 



xeAff 



xeA]^ 



which gives us 



fKM^)'dx< f K,i,{x)'dx<il + o{l))^ Yl 
•^0 Jo ^^^^ 



If we can prove that the function ^(dx, t) i-^ K^ijl(x, t^dxdt is continu- 

ous on D{[0, T],M), then E^{f^ K,ii{xfdx) converges to E°°{f^ K,fx{xfdx) 
for each N, and by the above inequality, this quantity is bounded by 



supjyE''^ Jq jj^xgAn -^t{x)'^dt , which is finite by hypothesis. Therefore 
the following two lemmas complete the proof of the proposition. □ 

Lemma 2.2. If lim sup E°° ( KfiJ,(xydxdt) < oo, then with probabil- 
ity one, IX is absolutely continuous, with ii{dx,t) — u{x,t)dx and 
E^{J^ J^u\x)dxdt) < oo. 
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Lemma 2.3. The function ij,{dx,t) K^fi{x,tydxdt is continuous on 



D{[0,T],M). 



Proof. First we prove Lemma 12.21 Let fk{x) = e^'^*'^^, and {fk}kez is an 
orthonormal basis for L^([0, 1]), with each continuous on T. For each 
measurable function g : [0,1] — )■ R, g{x)'^dxdt = \ {fk, g)\'^dt. 

We have (/o, K^fi) = 1 for all /i G Aii, and for k ^ 0, 



We conclude that as e goes to zero, | {fk, K^jj) \ increases to | (/fc, /i) |. Thus by 
the monotone convergence theorem. 



for all /i G D{\Q,T],M). Let iV = {/x G D{\Q,T],M,) : /o" Er=oo I (A, = 



oo}. Clearly Q°°{N) = by elementary measure-theoretical considera- 
tions. Thus, Y^'k'=oo\Uk, l^{t))\'^dt < oo Q°°-a..e. For such yU, the se- 
quence {fk,fi{t)) is in for almost all t, so there is a function u{x,t) = 
^"kL-ooifk^ such that Jq \\u\\ldt < oo. To see that jj{dx, t) = u{x, t)dx 

as a measure for each t such that u is finite, note that the collection {fk} 
forms an algebra which is dense in C(t). For each function g in the alge- 
bra, {g, tJi{t)) = {g, u{t)) and the measures must be the same. Now fi{dx, t) = 
u{x, t)dx for almost all t and therefore they are equal as functions in L'^{[0, T] x 
[0, 1]), as desired. 

To prove the continuity lemma, note that as — )■ z/ in A^i, K^fi{x) — >■ 
K^yi^x) when the endpoints are points of continuity of z/, that is, almost 
everywhere. Since K^^{xy is bounded by for /i G A^i, we can apply 
dominated convergence to see that Kf,n{x)'^dx — )■ Ky{x)'^dx. Then if 
/i. — )■ u. in D([0,T], A^), the time integrals converge, as desired. □ 




sin{2TTke) 



27rke 




17 



3 The weak form of the hydrodynamic equa- 
tion 

In this section we prove that the points of the hmiting measure satisfy a weak 
version of the hydrodynamic equation 10.21 First we identify the equation. 

Proposition 3.1. Let v be a solution to equations \0.2\0.(k Then the function 
u : [0, 1] X [0,T] R defined by 

u{x,t) = l{(_^^ty.x<s{t)}{v{x,t) + 1) 

satisfies the integral equation 



f{x)u{x,t)dx)— / f{x)u{x,0)dx 
Jo 

t rl ft pi 

f"{x)p{u{x,r))dxdr — a / / f'{x)u{x,r)dxdr (3.1) 
'o ^0 Jo Jo 

for each t <T. 

Proof. This is a calculation by integration by parts. By hypothesis, condi- 
tions [02E6] hold, and /'(O) = 0. Define v = ioi x > s{t). 

'(*) rm 

f{x)v{x,t)dx — / f{x)v{x,0)dx 




Jo 

t p^i^) 

dr / f{x)v{x,r)dxdr 
Jo 

t ft /•s{r) 

s'{r)f{s{r))v{s{r),r)+ / / f{x)drv{x,r))dxdr 
Jo Jo 

t ps{r) 

f{x){vr,r,{x,r) + avr,{x,r))dxdr 

Jo 

* s(r) r^''^ 

+ Of )|^i!q — a / f'{x)v{x,r)dxdr 
Jo 

s(r) r^'^ 

/'(x)f (x, r)|^i^Q + / f"{x)v{x.r)dxdr 
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' f{s{r)){-a~s'{r)) + af{0) 

+ / {f"{x) — af'{x))v{x,r)dxdt 
Jo 

s{t) r-s 

—af'{x)dx — dr j {r)f{x)dx 
Jo JO 

l-s{t) 

+ / {f"{x) - af{x))v{x,r)dxdt 

JO 

s(t) «(0) 

f{x)dx + / f{x)dx 
Jo 

+ / / f"{x)v{x,t)dx - / f\x){v{x,t) + l)lx<s{t)dxdt. 
Jo Jo Jo 

If it can be shown that f > in for fo > 0, foUowing our intuition for 
the heat equation, then the proposition is proven, since p{u) = v. To see 
that this holds, consider the space transformation T{x) = e°^, and define a 
function w in T{Dt) by 

w(T{x),t) = T\x)v{x,t). 

Then w satisfies Wt = a?y'^Wyy in T{Dt), and by Theorem 2.1 in |7], if w has 
a negative minimum, it occurs on the boundary, which must be at (l,to) foi" 
some to > because of the other boundary conditions. Then by Theorem 
2.14 in the same book, Vy{l,to) > 0. However, the boundary conditions 
require Wy(l,to) = —a, a contradiction. Therefore w > on T{Dt), and 
v>Oon (Dt). □ 

Next we turn to convergence. Our goal is the following lemma, leaving 
only uniqueness of such solutions to prove Theorem (0.1): 

Lemma 3.1. //, for each N, under are as in Lemma \2.1\ and the 
corresponding converge to the delta measure on an absolutely continuous 
measure with density uq G , then under Q'^ , for < t < T, fit{dx) is 
absolutely continuous a.s. with density u{x,t) such that, for f G C^([0,1]) 
such that /'(O) = /'(I) = 0, 

(/,«(-,t))-(/,«(-,0))= [\f",p{u{;s)))ds-a [\f',u{;s))ds. 
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Proof. The proof is by elementary estimates. We have aheady shown that for 
/ e C{[0, 1]), the P^-martingale Mt = {f,Xt)N - CN{f,Xs)Nds satisfies 

[|M,-Mo|>e]<-^, 

and 

CN{f,Xs)N = {ANf,p{Xs))N - a{DNf,Xs)N, 

at ([L5]) and respectively. For / G ^^([0, 1]) such that /'(O) = /'(I) = 0, 
Ajv/ — 7- /" and Dj^f — )■ /' uniformly, so 

/ {ANf,p{X,))N - a{DMf,Xs)Nds - [ {f",p{Xs))N-a{f',Xs)Nds 
Jo Jo 

goes to zero uniformly in D{[0,T], JUi). Now —a{f',Xs)N = —ci{f',Ps) in 
distribution, under and respectively, the right hand side being a 
continuous functional on A^i. However, the term with p is not one, so we 
must make some additional estimates. For p E Ai, let 



M,p{x) = ^ j ^\x-y\<edp{y). 



Because of the nature of the proof, this time we extend to the real line rather 
than the torus when integrating over y beyond the boundary points. In 
addition, extend / beyond [0, 1] so that /" is continuous and bounded, and 
interpret {f",p{M^pri)) to mean the integral over the whole real line (in this 
case, just [— e, 1 + e] because of the support of M^p). For the process X, 
recall that St = mm{x e An : Xt = 0}, defined to be (2A^ + 1)/2N for the 
state where all sites have one particle. Under P^, Xf G for all t a.s., 
so for X < St, Xt{x) > 1. Let Rt = St — ^/N, the last site with a particle. 
We claim that for < t < T, given 6 > 0, for e small enough, as goes to 
infinity. 



Q 



N 



{f,Pt)-{f.po)-{ / {f",piM,Ps))ds~a / {f,pMs 



> 6 



0. 



Because of the martingale bound and uniform convergence above, it suffices 
to show that for small enough e, 

i(r,p(M,/.,))-(r,^,(,,))i 
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vanishes as oo. Define 

h{x) = M,/i^(a;) - p(M,/i^(x)), 

and note that < h{x) < 1 and h{x) = for x > Rt + e. When Rt > 2e, for 
X & [e, Rt — e) , rjy > 1 for y G An such that ||/ — si < e. So for such x, 



2eN 



E 



> 



yeANn{x:\x-y\<e} 

\An n {a; : \x — y\ < e}\ 



2eN 



>- 



2eN-l 
2eN ' 



(3.2) 



and on this interval, e{x) := 1 — h{x) < l/(2eN). We will use the decompo- 
sitions 

f"{x)p{M,fir^{x))dx = j f"{x)M,fir^{x)dx - j ' f"{x)dx 



+ 



f"{x)e{x)dx + If 



where represents the left over bits of the integrals of bounded functions. 
Also, we have 

I f" {x)d^ip(^n.){x) = I f"{x)dfi^{x) - ^ ^ f"{x). 



Let M be such that |/"| < M, then by ([33), between e and - e, e{t) < 
l/2eN. Therefore we have 



Rt-e 



f"{x)e{x)dx 



< 



M 

2dV' 



As a Riemann sum, 



rRt 1 

/ nx)dx-- Yl /"( 



< c 



N, 
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for some going to zero, and < 4Me. Finally, in the space of distribu- 
tions, M^ji — )■ /i, which gives us convergence of (/", Mg/x) to (/", /i) uniformly 
on A^i, so we bound the difference by — 0. Combining estimates, we get 



i(r,p(M,/x,))-(r,/x,(,.))i < 



M 
2dV 



+ 4Me + Cn + a. 



Given 5 > 0, for all e small enough, as goes to infinity, the right hand 
side is less than and after integrating over t and taking expectations, the 
claim is proved. 

Next, J^{f" , p{M^fis))ds is a continuous functional on D{[0,T]Ai). The 
proof is the same as for in Lemma 2.3. Therefore, 



(/,/it)-(/,/io)- / {f",p{M,p,))ds~a / {f',fis)ds 



> 6 



0, 



for small e. Finally, p G -^^^([0, 1] x [0,T]) Q°°-a..s., and for such p, M^p jj, 
in L^, so 



{f",p{M,p,))ds- / {f",p{p,))ds -^0 
Jo 

a.s. Convergence is uniform on A^i fl L^, so 

{f,Pt) - {f,Po) - I {f",p{ps))ds- a {f',Ps)ds 
\Jo Jo J 

is less than 5 a.s. for all 5 > 0, and the lemma is proved. □ 



4 Uniqueness of weak solutions 

In the previous section, we found that a subsequential limit is concen- 
trated on solutions to the equation fl3.1l) . We need a uniqueness theorem for 
such u. 

Lemma 4.1. Given functions Uq G L'^{[0, I]), a G C([0,T]), a function u G 
L^([0, 1] X [0, T]) that satisfies the integral equation for all f G C^([0, 1]) 
with f'ifS) = /'(I) = is unique. 
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Proof. The proof is based on the method in A2.4 of [0], using techniques 
analagous to those used for the bound 12.11 Suppose u and u' are two 
solutions for a given uq, and let Ut{x) = u{x,t) — u'{x,t) and Pt{x) = 
p{u{x,t)) — p{u'{x,t)). Then 

dt{f,ut) = {f",p,)-a{f,%). (4.1) 

Again let ipo{x) = 1, ^pkix) = a/2 cos(7rA;a;) for A; > 0, and(f)k{x) = \/2sm{nkx). 
Then {ijjk}kLo is an orthonormal basis for L'^{[0, 1]) (the eigenf unctions for the 
Neumann problem). Note that {ipo.Ut) = for all t, and let 6fc(t) = {ipk.Ut). 
For positive integer iV, define 

fc=i 

a positive, differentiable function with 

dMt) = -2n' ^ hit) (^,, p,) - 27ra(t) ^ . 
fc=i fc=i 

Expanding in terms of {V'a;}, we get 

k ~2^2^ k 

k=l k=l j=l 

Define a{k,j) to be 1 when k — j is odd and otherwise, then 

2kcr{k,j) 

defined to be for k = j. Now we have 

k=l k=l j=l 

Absolute convergence of these sums to an L^([0,T]) function will allow us to 
use dominated convergence. By Schwarz's inequality, 

oo 

^\bk{s){^IJk,Ps}\ < \\Ut\\m[0,l])\\pt\\L2{[0,l]) < ll^t||i2([o,i]), 
k=l 
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which is in L^([0,T]) by hypothesis. Observe that - > P{k - jf, 
and for all k. 



< 



E4- 

Z — / 



m=l 



Then, 



EE 

k=i j=i 



hk{s)hj{s)a{k,j) 



oo oo 

EWE 

fc=l j=0 



bjcr{k,j) 



^2 _ j2 

sEi^'ilE'l 



k=l 

oo 



0=0 

oo \ 1/2 



0=0 



1/2 



fc=l 



0=0 



J^2 Z-^ jji2 
m=l 



C\\ut\\LH[o,i])J2j 



k=l 



<C\\ 

We apply dominated convergence to conclude that 



fc=l 



^2 _ j2 

27r2^6fc(s)(^fc,p,)rfs. (4.2) 



fc=i j=i 



A:=l 



Now, 



k=l 

since p is increasing. Since -R(O) = by hypothesis, R(t) < 0, so R(t) = 
for all t > 0, each bk(t) = 0, and Ut = 0, as desired. □ 

We can now prove Theorem 10. 1[ 

Theorem 4.1. Suppose that for each , Xq is a random vector with val- 
ues a.s. in Q^j- such that sup^ -E^[l/A^^^g^^ X^(a;)^] < oo, and that its 
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empirical measures {/i^((ix)}iv converge weakly to the delta measure on a 
fixed absolutely continuous measure UQ{x)dx. Then a limiting measure Q°° 
of the corresponding measures on D([0,T], A^) exists and is the delta 
measure on the unique solution of a weak version liS. 1\) of the Stefan problem 
( IO.^)) -/ f07^) with initial data uq. If a solution for the problem l\0.^) - l[07d\) with 
initial data {vq,so) = (p(Mo),inf{x : Uo{x) = 0}) exists, then the empirical 
measures of the process Yt converges to that solution in probability. 

Proof. The first part of tlie tlieorem follows from Lemmas ll.2[ I2.H 13. ![ and 
14.11 Given a strong solution (f , s) with initial data (p(mo), inf{a; : Uq{x) = 
0}), we have ^0(2;) = p{uq{x)) + l[o,so], which is easily checked since e 
a.s. Therefore, by Proposition 13. u{x) := v{x) + l[o,s(t)] is a solution 
to (13. ip and, by uniqueness, is the density of the limiting measure of the 
process Xf. In order for Yt = p{Xt) to converge to v = p{u), we need 
Pp{Xt) p{u{x,t))dx. Indeed, this fact is proved in the proof of Lemma 3.2, 
and the second part of the theorem is proved. □ 

5 Elastic exclusion 

Finally, we return to the regime of reflecting intervals by describing the hy- 
drodynamic limit of the exclusion process Zf. In this section we presume 
existence of a strong solution, a reasonable assumption that makes the fol- 
lowing proofs and calculations more natural. First we describe the dynamics 
of Zt more precisely. 

Each site is occupied by zero or one particles, except for 1 — 1/2N, which 
is always empty. Given Zt = 6, a state ( is accessible in one of three cases: 
first, if for some 1 < k < N - 1, d(2k-i)/2N = l,0{2k+i)/2N = 0, C(2k-i)/2N = 
0, C(2fc+i)/2Af = 1, and otherwise the states are equal. Second, if for A; = A^ — 1, 
we instead have C(2A;+i)/2Ar = 0, so that the particle is "killed" here. In either 
of these case, 6 moves to ( with rate iV^ times the length of the block of 
occupied sites to the left of {2k+l)/2N in 9. Third, if x is the first unoccupied 
site of 9, the state ( such that Cx = 1 and (y = 9y elsewhere is reached at 
rate aA^, and we think of the entire block being pushed over to make room 
for a new particle. For each A^, define, for x G An,x < St, 

Utix)=x+j^J2p(^^('^^- 

z<x 
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Note that Ut{St -1/N) = 1- 1/2N. Define 



OiiyeUtiA^n[0,St)) 
1 else 



so that, for x < St, between unoccupied sites y = Ut{x — 1/N) (or y = if 
X = 1/2N) and y' = Ut{x) of there are p{Xt{x)) occupied sites. It is 

not hard to check that the dynamics of are identical to those of Zt. 

We can write the inverse map explicitly: For y & An such that Zt{y) = 0, 
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Ttiy)^Ui\y) = y--}_^Ztiz) 

z<y 



and, given Zf, we extend the domain of Tj to all of with the same formula. 
We also define analagous transformations for the solution v of f l0.2p -f l(rB]) . Let 

Vt{x) = X + / v{x,t)dx, 
Jo 

for < X < s(t), and, since v is nonnegative and differentiable on [0,s(t)), 
we can define a differentiable inverse Tj = v^'^, and let 

z{y,t) ■■= ilj{v){y,t) := v{Tt{y),t)Tt{y), 

so that 

ry 

n{y) = y- z{y,t)dy. 
Jo 

and 

v{x,t) = z{vt{x),t)v[{x). 

One can check that for z so defined, < z < 1 on [0, 1]. Given z : [0, 1] — >■ 
[0, 1), we can similarly define r and v with the last two equations, where v is 
defined to be the inverse to r. Next we identify the hydrodynamic equation 
that functions z should satisfy: 

Lemma 5.1. Given functions z and v, defined on the appropriate regions, 
satisfying the relations above for all t, z is a solution to the differential equa- 
tion: 

dtz{y, t) = d,{K{z{y, t))dyz{y, t)) < y < 1, t > 0, (5.1) 

z{x,0) = zo{x) 0<x<l (5.2) 

K{z{y, t))dyz{y, t) \y=^ = a t > 0, (5.3) 

^(0,t) = t>0, (5.4) 
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with K{z) = 1/(1 — z)"^ , if and only if v is a solution to 1 0. ^If076l with s{t) 
1 — Jq z{x, t)dx. 

Proof. We will expand the equation 

dxxv{x, t) + adxv{x, t) = dtv{x, t), 

using 

dxVt{x) = 1 + v{x,t) 



1 - z{vt{x),t)' 

Let z = z{vt{x),t), let Zt and Zy refer to the partial derivatives of z with 
respect to its variables evaluated at {vt{x),t). We first check the equivalence 
of the conditions at the left boundary. The function v is continuous up to 
the boundary so v has a derivative at and the above equation holds. This 
gives us equivalence of 

Vx + a{y + 1) = 0, 
—z,, a 



(1-^)3 l-Z 



-Zy 



{i-zY 

aX. X = vt{x) = 0. Next, for < x < s(t), < t; < 1, 



Vx{x,t) - ^^3' 



Zy -1 



(1 - Z)2 \-Z 

-dx ,^ + 



l-z \l-zY {l-zY 



so 



Vxx + 



{l-z} 



vt{x) 



+ 



az„ 



(1-Z)2 (1-^)5 (I-Zf 
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Then 



dtVt{x) = - dtv{r,t)dr 
Jo 

= Vx{x, t) + av{x, t) — (fa.(0, t) + Of (0, t)). 
Under either set of conditions, the last term is equal to —a, giving 

dtVt{x) 



Z',1 



so 



dtv{x,t) 



(1-Z)3 

-dtVt{x)zy - zt 



1-z' 



+ 



az„ 



'yl-zf ' {l-Z^ (1-^)3' 

and Vxx + = Vt is equivalent, for y = Vt{x), to 

dyz{y,t) 



dt{z{y,t)) = dy 



z{y,t)y 



as desired. 



□ 



Theorem 5.1. The process Zt = \l/(Xt), where Xt has initial data uq as in 
Theorem 4jJ_. converges weakly in D{[0,T], M.) to the measure with density 
that is the unique solution to the differential equation with initial 

data Zo = iP{uq), assuming such a solution exists. 

Proof. We will prove that for 5 > 0, / G C([0, 1]), 



pN 


sup 


L 




.0<t<T 





f{y)f^Zt(y)dy- / f{y)z{y,t)dy 



> 6 



->0, 



which suffices to prove the theorem. First, since v is sufficiently smooth, if 
we extend f{vt{x)) to be where it is not defined, then our previous theorem 
gives 



pN 


sup 






.0<i<T 





f{vt{x))i2p(^x)t{x)dx- / f{vt{x))v{x,t)dx 



> 6 
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and since the last terms in the two differences are equivalent by a change of 
variables, we need to look at the difference of 



1 

/ fiy)f^zt(y) = T7 f{y)Zt{y) 



and 



° xgAjv y€AN 

Since Z < 1 and / is continuous on a compact interval, the problem is 
reduced to the difference of y and VtiTtiy)). By Theorem 14.11 and the fact 
that the limiting measure is continuous, we have 



sup sup 

o<t<ro<a<i 



X 



> 6 



[E] ^ 0, 



where E is the set on the left hand side. For fixed t, suppose y = Vt{x) = 
X + Jq v{z,t)dz, and we claim that for large enough A^, in the set E'^, 

\Tt{y)-x\<6. 

Indeed, suppose Tt{y) < x — 6. Then 

y-j^Yl <x-6 

z<y 

^ J2 PiXt) > I v{z,t)dz + 5 

z<Tt{y) 

z<Ttiy) z<x 

and Tt{y) > x, a contradiction. A contradiction results from the opposite 
inequality in the same way. Then, vt being continuous, we obtain a bound 
on the difference 

\MTtiy)) - vt{x)\ = \vt{Tt{y)) - y\ , 
which gives the result we need and completes the proof of the theorem. 

□ 
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